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Abstract 

Of concern in this paper is an investigation of heat and mass transfer over a stretching 
sheet under the influence of an appUed uniform magnetic field and the effects of Hall cur- 
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J>-,' rent are taken into account. The non-linear boundary layer equations together with the 

^ ■ 
O ' boundary conditions are reduced to a system of non-linear ordinary differential equations 

by using the similarity transformation. The system of non-linear ordinary differential 

^ . equations are solved by developing a suitable numerical techniques such as finite differ- 



ence scheme and Newton's method of linearization. The numerical results concerned with 
the velocity, temperature and concentration profiles as well as the skin-friction coefficient, 
local Nusselt number Nu and the local sherhood number Sh for various values of the non- 
dimensional parameters presented graphically. 

Keywords: Hall current, Magnetohydrodynamic, Stretching sheet, Skin-friction, Heat 
transfer rate. Chemical reaction 

1 Introduction 

The study of MHD free-convective flow and mass transfer over a stretching sheet has become a 
considerable attention to the researchers due to its many engineering and industrial applications 
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such as, in polymer processing, electro-chemistry, MHD power generators as well as in flight 
magnetohydrodynamics. In the extrusion of a polymer sheet from a die, the sheet is some times 
stretched. By drawing such a sheet in a viscous fluid, the rate of cooling can be controlled 
and the final product of the desired characteristics can be achieved. This problem has also an 
important bearing on metallurgy where magnetohydrodynamic (MHD) techniques have recently 
been used. 

Crane (Crane 1970) first introduced the study of steady two-dimensional boundary layer 
fiow caused by a stretching sheet whose velocity varies linearly with the distance from a fixed 
point in the sheet. Later on several investigators (cf. Gupta and Gupta 1977; Rajagopal et al. 
1984; Siddapa and Abel 1985; Chen and Char 1988; Laha et al. 1989; Vajravelu and Nayfeh 
1992; South et al. 2002 and Tan et al. 2008) studied various aspects of this problem, such as 
the heat, mass and momentum transfer in viscous fiows with or without suction or blowing. 
The infiuence of a uniform transverse magnetic field on the motion of an electrically conducting 
fiuid past a stretching sheet was investigated by Pavlov (Pavlov 1974), Chakraborty and Gupta 
(Chakraborty and Gupta 1979), Kumari et al. (Kumari et al. 1990), Andersson (Andersson 
1992, Andersson et al. (Andersson et al. 1992) and Char (Char 1994). The effect of chemical 
reaction on free-convective flow and mass transfer of a viscous, incompressible and electrically 
conducting fluid over a stretching sheet was investigated by Aflfy (Aflfy 2004) in the presence of 
transverse magnetic fleld. In all these investigations the electrical conductivity of the fluid was 
assumed to be uniform. However, in an ionized fluid where the density is low and/or magnetic 
fleld is very strong, the conductivity normal to the magnetic fleld is reduced due to the spiraling 
of electrons and ions about the magnetic lines of force before collisions take place and a current 
induced in a direction normal to both the electric and magnetic flelds. This phenomenon avail- 
able in the literature, is known as Hall effect. Thus, the study of magnetohydrodynamic viscous 
flows, heat and mass transfer with Hall currents has important bearing in the engineering ap- 
plications. 

Hall effects on MHD boundary layer flow over a continuous semi-inflnite flat plate moving 
with a uniform velocity in its own plane in an incompressible viscous and electrically conducting 
fluid in the presence of a uniform transverse magnetic fleld were investigated by Watanabe and 



Pop (Watanabe and Pop 1995). Aboeldahab and Elbarbary (Aboeldahab and Elbarbary 2001) 
studied the Hall current effects on MHD free-convection flow past a semi-infinite vertical plate 
with mass transfer. The effect of Hall current on the steady magnetohydrodynamics flow of an 
electrically conducting, incompressible Burger's fluid between two parallel electrically insulating 
infinite planes was studied by Rana et al. (Rana et al. 2008). 

The aim of this paper is to study the Hall effects on the steady MHD free-convective flow 
and mass transfer over a stretching sheet in the presence of a uniform transverse magnetic field. 
The boundary layer equations are transformed by a similarity transformation into a system of 
non-linear ordinary differential equations and which are solved numerically by using the finite 
difference technique. Numerical calculations were performed for various values of the magnetic 
parameter, Hall parameter and the relative effect of chemical diffusion on thermal diffusion pa- 
rameters. The results are given for the velocity distribution and the coefficients of skin-friction 
along and transverse to the direction of motion of the stretching surface. Similarity solution for 
the temperature field in the above flow is also found and the rate of heat transfer at the stretch- 
ing sheet is computed for various values of magnetic and Hall parameters. Such a study is also 
applicable to the elongation to the bubbles and in bioengineering where the flexible surfaces of 
the biological conduits, cells and membranes in living systems are typically lined or surrounded 
with fluids which are electrically conducting (e.g., blood) and being stretched constantly. 



2 Mathematical Formulation 

We consider the steady free-convective flow and mass transfer of an incompressible, viscous and 
electrically conducting fluid past a flat surface which is assuming from a horizontal slit on a 
vertical surface and is stretched with a velocity proportional to distance from a fixed origin O 
(cf. Fig. 1). We choose a stationary frame of reference {x,y,z) such that x-axis is along the 
direction of motion of the stretching surface, y-axis is normal to this surface and z-axis is trans- 
verse to the xy-plane. A uniform magnetic field Bq is imposed along y-axis and the effect of Hall 
currents is taken into account. The temperature and the species concentration are maintained 



at a prescribed constant values T^, C^ at the sheet and Too and Coo are the fixed values far 
away from the sheet. 

Taking Hall effects into account the generalized Ohm's law (cf. Cowling (Cowling 1957)) 
may be put in the form: 



J = -^-^ (e + VxB- —J X b] , 
1 + m"^ V eue / 



in which V represents the velocity vector, E is the intensity vector of the electric field, B is 
the magnetic induction vector, H the intensity of the magnetic field vector, /Xg the magnetic 
permeability, J the electric current density vector, m = ^^ is the Hall parameter, a the elec- 
trical conductivity, e the charge of the electron and Ug is the number density of the electron. 
The effect of Hall current gives rise to a force in the ^-direction which in turn produces a cross 
flow velocity in this direction and thus the flow becomes three-dimensional. To, simplify the 
analysis, we assume that the flow quantities do not vary along z-direction and this will be valid 
if the surface is of very large width along the z-direction. 

Under these assumptions the boundary layer free-convection flow with mass transfer and 
generalized Ohm's law is governed by the following system of equations: 

dx dy 



du du d'^u „ ,^ ^ , n*/^ ^ \ (^Bn . , .^. 



dw dw d^w gBI . 

dx dy dy"^ p(l -|- m?) 



dT dT d'^T 



dx dy dy"^ ' 



where [u^v^vS) are the velocity components along the {x,y,z) directions respectively, u is 
the kinematic viscosity, g the acceleration due to gravity, /3 the coefficient of thermal expansion, 
P* the coefficient of expansion with concentration, T and C are the temperature and concen- 
tration respectively, p the density of fluid, a the thermal diffusivity, D the thermal molecular 
diffusivity, kQ the reaction rate constant and n the order of reaction. 

The boundary conditions for this problem can be written as 

u = bx, V = w = 0, T = T^, C = Cw at y = (6) 

u = w = 0, T = Too, C = Coo as y ^ oo (7) 

where b > 0. The boundary conditions on the velocity in (6) are the no-slip conditions at the 
surface y = 0, while the boundary conditions on velocity at y— )■ oo follow from the fact that 
there is no flow far way from the stretching surface. 

The temperature and species concentration are maintained at a prescribed constant values T^ 
and Cw at the sheet and are assume to vanish far way from the sheet. 

Introducing the similarity transformation: 

u = bxf'(r]), V = —vbufirj), w = bxgir]),-)] = \ —y (8) 

V ly 
T — T r — r 



T —T n — c 

where a prime denotes derivative with respect to rj. Substituting (8) into equations (2)- (7), 
yields 

/'" + //" - r + GrO + Gc<t> - -^ (/' + mg) = (9) 

1 + m^ 



/ + /*'-(/' + T^)<; + T^/' = o (10) 

V 1 + m^ / 1 + m^ 

e" + Prfe' = (11) 

0" + 5c(07-70") = O (12) 

where M = ^-^ is the magnetic parameter, 
Qrf, _ gp( -w- oo) ^]-^g local Grashof number, 

Gc = ,1?~ °° the local modified Grashof number, 

m = ^^ is the Hall current parameter, 

Pr = - the Prandtl number, 

Sc= jj the Schmidt number, 

and 7 = ^(Ct„ — Coo)"^^ the non-dimensional chemical reaction parameter. 

The boundary conditions (6) and (7) are now obtained from (8) as 

/'(0) = 1, /(O) = (?(0) = 0, ^(O) = 0(O) = 1 (13) 

/'(oo) = gioo) = ^(oo) = 0(oo) = 0. (14) 

The major physical quantities of interest are the skin- friction coefficient Cf, local Nusselt number 
Nu and the local Sherwood number Sh are defined respectively by, 

Cf = -^ = no), where r^ = ^ (^] = ^.[^axf"{0), (15) 

MA/fax \9y)y=Q V/y 



Nu "" ^"-' '-- - ''"^^^ 



^'(0), where q^ = -k{J-\ = -k^(T^ - T^)9\0), (16) 



kMT^-T^) \dy) 



Sh '" 



-(/)'(0), where m^ = -D l^\ = -D^{C^ - Coo)0'(O), (17) 
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When M = m = and Gr = Gc = 0, present flow problem becomes hydro dynamic 
boundary layer flow past a stretching sheet whose analytical solution put forwarded by Crane 
(Crane 1970) as follows: 



f{r^) = l-e'\ I.e. f'{v) 



e 



-V 



An attempt has been made to validate our results for the axial velocity f'{rj), we compared our 
results with this analytical solution. 



3 Numerical Methods 

Several authors Andersson et al. (Andersson et al. 1992), Aflfy (Aflfy 2004) and Abo-Eldahab 
(Abo-Eldahab 2001) used numerical techniques for the solution of such two-point boundary 
value problems is the Runge-Kutta integration scheme along with the shooting method. Al- 
though this method provides satisfactory results, it may fail when applied to problems in which 
the differential equations are very sensitive to the choice of its missing initial conditions. More- 
over, difficulty arises in the case in which one end of the range of integration is at infinity. The 
end point of integration is usually approximated by replacing a finite representation to this point 
and it is obtained by estimating a value at which the solution will reach its asymptotic state. 
On the contrary to the above mentioned numerical method, we used in the present paper has 
better stability, simple, accurate and more efficient numerical technique. The essential features 
of this technique is that it is based on a finite difference scheme with central differencing and 
based on the iterative procedure. 

We substitute F = f and equations (9) and (10) are then transformed as 

F" + fF' -F^ + Gr9 + Gc(f) ^^ (F + mg) = (18) 

1 + rri'^ 

»" + /»'- (f + T^)9+^F = (19) 

V 1 + m^ / 1 + m"^ 

with the boundary conditions 

F(0) = 1, F(oo) = 0, /(0)=0, (7(0) = 0, g{oo) = (20) 
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Central difference scheme is used for derivatives with respect to rj : 



Wi+i - Wi-x , ^//^„^2^ 



(«'.). = "2^^ + ^((^^^ 



where w stands for F, g, 6 and 0; i is the grid index in the r7-direction with 



i'^^v)^ = 1-2 + OiiSrj) 



5rj being the increment along 77-direction. 

We use Newton's linearization method to linearize the discretised equations as follows. We 
assume that the values of the dependent variables at the /cth iteration are known. Then the 
values of the variables at the next iteration are obtained from the following equation 

«,f+i = u,f + 5wf (21) 

where 5w^ represents the error at the k-th iteration. Using (21) in (18) nd (19) and drop- 
ping quadratic terms in ^tof , we get a system of linear algebraic equations for (5t«f . The re- 
sulting system of block tri-diagonal equations is then solved by Thomas algorithm. With an 
aim to test the accuracy of this numerical method we have compared the values of f'{rj) for 
M = m = Gr = Gc = with those of analytical studies of Crane (Crane 1970) and have found 
excellent agreement presented through the Fig. 2. 



4 Results and Discussion 

The system of ordinary differential equations (9)-(12) subject to the boundary conditions (13) 
and (14) are solved numerically by employing a finite difference scheme with Newton's lineariza- 
tion method described in the previous section. For numerical computations the following values 
of the physical parameters have been considered according to the data used in (Afify 2004): 

Pr = 0.7, Gr = 0.5, Gc = 0.5, Sc = 0.5, n = 1, 2, 3; M = to 5.0; m = to 3; 



7 = 0.1, 0.5, 1.0 

Figs. 3-5 depicted the variation of velocity profiles /'(?]), frj) and g{r]) corresponds to the ax- 
ial velocity, transverse velocity and the velocity along z-direction called cross fiow velocity for 
different values of M. From these three figures we see that the velocity decreases with the 
increase in the magnetic parameter M. This is to be expected from the physical consideration 
since as M increases, Lorentz force which opposes the fiow and leads to deceleration of the fiuid 
motion. It is interesting to note that for a fixed values of M and m, /{rj) reaches a uniform 
value asymptotically at a certain height t] above the sheet. By contrast, the cross fiow velocity 
component induced due to Hall effects and shows a anomalous behaviour with variation of M. 
Fig. 5 shows that for fixed values of M and m, g{r]) reaches a maximum value at a certain 
height 1] beyond which it decreases gradually in asymptotic nature. It can be noted from Fig. 5 
that in the absence of magnetic parameter [M = 0), cross fiow velocity g{7]) vanishes. 

Figs. 6 and 7 illustrated distribution of non-dimensional temperature 6{ri) and concentra- 
tion (j){rj) for different values of magnetic parameter M with a fixed values m = 0.1, Pr = 0.7, 
Gr = Gc = 0.5. These two figures shows that the temperature as well as concentration increases 
with with the increase of magnetic parameter M. It may conclude that in the presence of mag- 
netic field temperature and species concentration can be increased. 

Effects of Hall current parameter m on the velocity profiles are presented in Figs. 8-10. For 
a fixed value of M the velocity components /'(?]) and f{r]) at a given height increases monoton- 
ically with the increase in m, while the trend is reverse in the case of cross-fiow velocity g{ri) 
induced by Hall effects. It is interesting to observe from Figs. 8 and 9 that velocity greatly 
affected for m < 3, beyond which no significant change occur. Figs. 11 and 12 noticed that 
the temperature and concentration profiles also affected by the Hall current parameter m. It is 
shown from Figs. 11 and 12 that the temperature as well as concentration decreases with the 
increase in Hall parameter m. 

The effects of chemical reaction on free-convective fiow and mass transfer of an electrically 
conducting fiuid over a stretching sheet have also been studied in the presence of magnetic field. 
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Figs. 13 and 14 shows that velocity and temperature shghtly decreases with the increase of 
the chemical reaction parameter 7, whereas, the concentration significantly diminishes with an 
increase in 7. 

Fig. 16 shows the variation of skin-friction with M for different values of Hall parameter m. 
The skin-friction increases with the increase of both magnetic parameter M and Hall parameter 
m. The rate of heat transfer at the sheet is increases with an increase in m but reduces with 
increasing M shown in Fig. 17. Variation of the rate of mass transfer at the wall is presented 
in Fig. 18, which shows that the mass transfer rate decreases with M but increases with the 
increase of m. It is interesting to note that the skin-friction coefficient, rate of heat transfer and 
mass transfer at the sheet vary linearly with M. 



5 Conclusions 

The effects of Hall current parameter on free-convective flow and mass transfer of a viscous, 
incompressible and electrically conducting fluid over a stretching sheet have been studied in the 
presence of magnetic fleld. The non-linear boundary value problem is solved numerically by 
flnite difference scheme with Newton's linearization method. 

From the present investigation, it may be concluded that all the instantaneous flow charac- 
teristics are affected by the Hall current parameter m. In the presence of an external magnetic 
fleld of an electrically conducting fluid. Hall current induced which in turn produces a cross flow 
velocity in the direction perpendicular to the directions of both axial and transverse velocity. 
The velocity proflles /', /, g decreases with the increase in M whereas the temperature as well 
as concentration increases with the increase of magnetic fleld strength. The flow velocities /' 
and / gradually increases with the increase of Hall current parameter m while, in the case of 
cross flow velocity g and the temperature as well as the concentration a reversal trend is observed. 

Acknowledgement: The author G. C. Shit is thankful to the UGC, New Delhi for the 
financial support during this investigation. 
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Fig. 1 Physical sketch of the problem 
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Fig. 2 Variation of f'{i]) with t] ioT M = m = and Gr = Gc = 
(Comparison of the present study with the analytical solution of Crane (1970) 
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Fig. 3 Variation of f'{i]) with t] for different values of M and 




Fig. 4 Variation of f{r]) with 77 for different values of M and 
Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, n = 1, m = 1, 7 = 0.1 
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Fig. 5 Variation of g{ri) with 7] for different values of M and 




Fig. 6 Distribution of non-dimensional temperature 6{rj) for different M with 
Pr = 0.7, Gr = Gc= 0.5, Sc = 0.5, n=l, m = l, 7 = 0.1 
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Fig. 7 The concentration profile (f){ri) for different M with 




Fig. 8 Effect of Hall parameter m on axial velocity f'{r]) for 
M = 1.0, Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, 7 = 0.1, n = 1 
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Fig. 9 Effect of Hall parameter m on axial velocity /{r]) for 

M = 1.0, Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, 7 = 0.1, n = 1 
0.12 




Fig. 10 Effect of Hall parameter m on the z-direction velocity §{7]) for 
M = 1.0, Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, 7 = 0.1, n = 1 
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Fig. 11 Effect of Hall parameter m on non-dimensional temperature ^(77) for 




Fig. 12 Effect of Hall parameter m on the concentration profile ^(77) for 
M = 1.0, Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, 7 = 0.1, n = 1 




Fig. 13 Variation of axial velocity profile /'(?]) with rj for different 7 and 

Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, M = 0.5, m = 0.1, n = 1 
1 
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Fig. 14 Variation of the temperature profile ^(77) with rj for different 7 and 
Pr = 0.7, Gr = Gc = 0.5, 5c = 0.5, M = 0.5, m = 0.1, n = 1 
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Fig. 15 Variation of the concentration profile 0(?7) for different 7 and 

Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, M = 0.5, m = 0.1, n = 1 
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Fig. 16 Variation of skin-friction for different values of Hall parameter m and 
Pr = 0.7, Gr = Gc = 0.5, ^c = 0.5, n = 1, 7 = 0.1 
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Fig. 17 Distribution of heat transfer rate at the channel wall for different values of Hall 

parameter m and Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, n = 1, 7 = 0.1 
0.49 




Fig. 18 Distribution of the rate of concentration at the channel wall for different values of Hall 
parameter m and Pr = 0.7, Gr = Gc = 0.5, Sc = 0.5, n = 1, 7 = 0.1 
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